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Abstract 



o 

-^ Polynomial algebra offers a standard approach to handle several problems in geometric mod- 

eling. A key tool is the discriminant of a univariate polynomial, or of a well-constrained system 
of polynomial equations, which expresses the existence of a multiple root. We describe discrimi- 
nants in a general context, and focus on exploiting the sparseness of polynomials via the theory 
of Newton polytopes and sparse (or toric) elimination. We concentrate on bivariate polynomials 
and establish an original formula that relates the mixed discriminant of two bivariate Laurent 
polynomials with fixed support, with the sparse resultant of these polynomials and their toric 

CT\ ■ Jacobian. This allows us to obtain a new proof for the bidegree of the mixed discriminant as 

well as to establish multipicativity formulas arising when one polynomial can be factored. 

in 

q ; 1 Introduction 

m 

Polynomial algebra offers a standard and powerful approach to handle several problems in geo- 
metric modeling. In particular, the study and solution of systems of polynomial equations has 
been a major topic. Discriminants provide a key tool when examining well-constrained systems, 
$_i . including the case of one univariate polynomial. Their theoretical study is a thriving and fruitful 



domain today, but they are also very useful in a variety of applications. 

The best studied discriminant is probably known since high school, where one studies the 
discriminant of a quadratic polynomial f(x) = ax 2 + bx + c = (a ^ 0). The polynomial / has 
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a double root if and only if its discriminant A2 = b 2 — Aac is equal to zero. Equivalently, this 
can be denned as the condition for f(x) and its derivative f'(x) to have a common root: 

3 x : f{x) = ax 2 + bx + c = f'{x) = 2ax + b = 44> A 2 = 0. (1) 

One can similarly consider the discriminant of a univariate polynomial of any degree. If we wish 
to calculate the discriminant As(/) of a polynomial / of degree five in one variable, we consider 
the condition that both / and its derivative vanish: 

f(x) = ax 5 + bx 4 + ex 3 + dx 2 + ex + g = 0, 
f'(x) = 5ax 4 + Abx 3 + 3cx 2 + 2dx + e = 0. 

In this case, elimination theory reduces the computation of A5 to the computation of a 9 x 9 
Sylvester determinant, which equals a As(/). If we develop this determinant, we find out that 
the number monomials in the discriminant increases rapidly with the input degree: 

A 5 = -2050a 2 g 2 bedc + 356abed 2 c 2 g - 80b 3 ed 2 cg + I8dc 3 b 2 g 
e - 746agdcb 2 e 2 + 144a6 2 e 4 c - 6ab 2 e 3 d 2 - 192a 2 be 4 d - Ad 2 ac 
3 e 2 + U4d 2 a 2 ce 3 - Ad 3 b 3 e 2 - Ac 3 e 3 b 2 - 80abe 3 dc 2 + 18b 3 e 3 
dc + 18d 3 acbe 2 + d 2 c 2 b 2 e 2 - 276 4 e 4 - 128a 2 e 4 c 2 + 16ac 4 e 3 - 27 
a 2 d 4 e 2 + 256a 3 e 5 + 3125aV + 160a 2 gbe 3 c + 560a 2 gdc 2 e 2 + 1020 
a 2 gbd 2 e 2 + 160ag 2 b 3 ed + 560ag 2 d 2 cb 2 + W20ag 2 b 2 c 2 e - 192 
b 4 ecg 2 + 2Aab 2 ed 3 g + 2Aabe 2 c 3 g + UAb 4 e 2 dg - 6b 3 e 2 c 2 g + 14 
4dc 2 &V - 630dac 3 bg 2 - 630d 3 a 2 ceg - 72d 4 acbg - 72dac 4 e 
g - Ad 3 c 2 b 2 g - 1600a3 3 c6 3 - 2500a 3 5 3 6e - 50a 2 5 2 6 2 e 2 - 3750a 3 
g 3 dc + 2000a 2 g 3 db 2 + 2000a 3 g 2 ce 2 + 825a 2 g 2 d 2 c 2 + 2250a V& 
c 2 + 2250a 3 c/ 2 ed 2 - 900aV&d 3 - 900a 2 ^ 2 c 3 e - 36agb 3 e 3 - 1600 
a 3 ge 3 d + 16d 3 ac 3 g - 128d 2 b 4 g 2 + 16d 4 b 3 g - 27c 4 b 2 g 2 + 108ac 5 
g 2 + 108a 2 d 5 g + 256b 5 g 3 . 

In fact, if we compute the resultant of / and xf by means of the 10 x 10 Sylvester determinant, 
we find the more symmetric output: agA^(f). This formula is very well known for univariate 
discriminants [T7], and we generalize it in Theorem [2j 

One univariate polynomial is the smallest well-constrained system. We are concerned with 
multivariate systems of sparse polynomials, in other words, polynomials with fixed support, 
or set of nonzero terms. Sparse (or toric) elimination theory concerns the study of resultants 
and discriminants associated with toric varieties. This theory has its origin in the work of 
Gel'fand, Kapranov and Zelevinsky on multivariate hypergeometric functions. Discriminants 
arise as singularities of such functions [18] . 

Gel'fand, Kapranov and Zelevinsky [17J established a general definition of sparse discrimi- 
nant, which gives as special case the following definition of (sparse) mixed discriminant (see Sec- 
tion [2] for the relation with the discriminant of the associated Cayley matrix and with the notion 
of mixed discriminant in [3]). In case n = 2, the mixed discriminant detects tangencies between 



families of curves with fixed supports. In general, the mixed discriminant A-Ai,...,A n (fli ■ ■ ■ > fn) 
of n polynomials in n variables with fixed supports A\, . . . , A n C 7L n is the irreducible polynomial 
(with integer coprime coefficients, defined up to sign) in the coefficients of the /, which vanishes 
whenever the system f\ = • • • = f n = has a multiple root (that is, a root which is not simple) 
with non-zero coordinates, in case this discriminantal variety is a hypersurface (and equal to 
the constant 1 otherwise). The zero locus of the mixed discriminant is the variety of ill-posed 
systems [23]. We shall work with the polynomial defining the discriminant cycle (see Section [2]) 

which is defined as the power A^ 1 '"7 of the mixed discriminant raised to the index 

i(A x , . ..,A n ) = [Z n : ZAi + • • • + Z,A n ], (2) 

which stands for the index of lattice ZA\ + • • • + 7jA n in Z n . In general, this index equals 1 and 
so both concepts coincide. 

Discriminants have many applications. Besides the classical application in the realm of 
differential equations to describe singularities, discriminants occur for instance in the description 
of the topology of real algebraic plane curves [19J, in solving systems of polynomial inequalities 
and zero-dimensional systems [16], in determining the number of real roots of square systems 
of sparse polynomials [8], in studying the stability of numerical solving [6], in the computation 
of the Voronoi diagram of curved objects [13], or in the determination of cusp points of parallel 
manipulators [20] . 

Computing (mixed) discriminants is a (difficult) elimination problem. In principle, they can 
be computed with Grobner bases, but this is very inefficient in general since these polynomials 
have a rich combinatorial structure [17] . Ad- hoc computations via complexes (i.e., via tailored 
homological algebra) are also possible, but they also turn out to be complicated. The tropical 
approach to compute discriminants was initiated in [10] and the tropicalization of mixed pla- 
nar discriminants was described in [9]. Recently, in [12] . the authors focus on computing the 
discriminant of a multivariate polynomial via interpolation, based on [111 123] : the latter essen- 
tially offers an algorithm for predicting the discriminant's Newton polytope, hence its nonzero 
terms. This yields a new output-sensitive algorithm which, however, remains to be juxtaposed 
in practice to earlier approaches. 

We mainly work in the case n = 2, where the results are more transparent and the basic 
ideas are already present, but all our results and methods can be generalized to any number of 
variables. This will be addressed in a subsequent paper [7]. Consider for instance a system of 
two polynomials in two variables and assume that, the first polynomial factors as fi = f[ ■ /{'• 
Then, the discriminant also factors and we thus obtain a multiplicativity formula for it, which 
we make precise in Corollary [6] This significantly simplifies the discriminant's computation and 
generalizes the formula in [2] for the classical homogeneous case. This multiplicativity formula 
is a consequence of our main result (Theorem [2] in dimension 2, see also Theorem [3] in any 
dimension) relating the mixed discriminant and the resultant of the given polynomials and their 
toric Jacobian (see Section [3] for precise definitions and statements). As another consequence 
of Theorem [2] we reprove, in Corollary [5] the bidegree formula for planar mixed discriminants 
in [3]. 



The rest of this chapter is organized as follows. The next section overviews relevant existing 
work and definitions. In Section [3] we present our main results relating the mixed discriminant 
with the sparse resultant of the two polynomials and their toric Jacobian. In Section[5]we deduce 
the general multiplicativity formula for the mixed discriminant when one polynomial factors. 

2 Previous work and notation 

In this section we give a general description of discriminants and some definitions and notations 
that we are going to use in the following sections. 

Given a set A C M n , let Q = conv(A) denote the convex hull of A. We say that A is a lattice 
set or configuration if it is contained in Z n , whereas a polytope with integer vertices is called a 
lattice polytope. We denote by Vol(-) the volume of a lattice polytope, normalized with respect 
to the lattice Z n , so that a primitive simplex has normalized volume equal to 1. Normalized 
volume is obtained by multiplying Euclidean volume by n\. 

Given a non-zero Laurent polynomial 



/ = £ 



CnX , 



the finite subset A of Z n of those exponents a for which c a ^ is called the support of /. The 
Newton polytope N(f) of / is the lattice polytope defined as the convex hull of A. 

A (finite) set A is said to be full, if it consists of all the lattice points in its convex hull. In [3], 
A is called dense in this case, but we prefer to reserve the word dense to refer to the classical 
homogeneous case. A subset F C A is called a face of A, denoted F -< A, if F is the intersection 
of A with a face of the polytope conv(A). 

As usual Qi + Q2 denotes the Minkowski sum of sets Q\ and Q2 in W 1 . The mixed volume 
MV(Q±, . . . ,Q n ) of n convex polytopes Qi in W 1 is the multilinear function with respect to 
Minkowski sum that generalizes the notion of volume in the sense that MV(Q, . . . ,Q) = Vol(Q), 
when all Qi equal a fixed convex polytope Q. 

The following key result is due to Bernstein and Kouchnirenko. The mixed volume of the 
Newton polytopes of n Laurent polynomials /i(x), ... ,f n (x) in n variables is an integer that 
bounds the number of isolated common solutions of fi(x) = 0, . . . , f n (x) = in the algebraic 
torus (K*) n , over an algebraically closed field K containing the coefficients. If the coefficients 
of the polynomials are generic, then the common solutions are isolated and their number equals 
the mixed volume. This bound generalized Bezout's classical bound to the sparse case: for 
homogeneous polynomials the mixed volume and Bezout's bound coincide. 

Mixed volume can be defined in terms of Minkowski sum volumes as follows. 

MV n (Q 1 , ...,Q n ) = ^(-l)"- fe Y, ^ Vol ( E Qi) ■ 

fe=l Ic{l,...,n},\I\=k ' iel 



This implies, for n = 2: 

2MV(Q 1 ,Q 2 ) = Vol(Qi + Q 2 ) - Vol(Qi) - Vol(Q 2 ). 

A family of finite lattice configurations Ai, . . . , A^ in Z™ is called essential if the affine di- 
mension of the lattice ZAi + • • • + TLA^ equals k — 1, and for all proper subsets I C {1, . . . , k} it 
holds that the affine dimension of the lattice generated by {Ai,i € 1} is greater or equal than 
its cardinality \I\. 

Definition/Theorem 1. [T7J[25] Fix a family ofn + 1 finite lattice configurations A±, . . . ,A n+ ± 
which contains a unique essential subfamily {Ai,i E /}. Given Laurent polynomials in n vari- 
ables /i, . . . ,/n+i with supports Ai, . . . , A n+ \, the resultant ResAi,...,A n+ i(fl, . . . , / n +i) is the 
irreducible polynomial with coprime integer coefficients (defined up to sign) in the coefficients of 
/i, . . . , fn+i, which vanishes whenever f\, . . . , / n +i have a common root in the torus (C*) n . In 
fact, in this case, the resultant only depends on the coefficients of fi with i £ I. 

If there exist more than one essential subfamilies, then the (closure of the) variety of solvable 
systems is not a hypersurface and in this case we set: 

R-eS J 4 1 ,..., J 4 n+1 (/l,...,/n+l) = 1. 

In what follows, we consider n (finite) lattice configurations A\, . . . , A n in Z n and we denote 
by Qi, . . . , Q n their respective convex hulls. Let /i, . . . , f n be Laurent polynomials with support 
Ai, . . . , A n respectively: 

fi{x) =S ^c i>a x a , i = l...,n. 

In [3 J the mixed discriminantal variety, is defined as closure of the locus of coefficients Cj iCf for 
which the associated system /i = • • • = f n = has a non-degenerate multiple root x G (K*) n . 
This means that x is an isolated root and the n gradient vectors 

dfj , > dfj 

dx\ ' ' dx n 

are linearly dependent, but any n — 1 of them are linearly independent. 

If the mixed discriminantal variety is a hypersurface, the mixed discriminant of the previous 
system is the unique up to sign irreducible polynomial AA 1) ... t A n with integer coefficients in the 
unknowns Cj jfl which defines this hypersurface. Otherwise, the family is said to be defective and 
we set &-Ai,...,A n = 1- The mixed discriminant cycle A.Ai,...,A n is equal to i(A±, . . . , A n ) times the 
mixed discriminant variety, and thus its equation equals /S.A x ,...,A n raised to this integer (defined 
in fll])). 



By O Theorem 2.1], when the family A\, . . . ,A n is non defective, the mixed discriminant 
A-Ai,...,A n coincides with the A-discriminant defined in |17j . where A is the Cayley matrix 
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( 1 








1 



Ai 



o \ 
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Am 



J 



This matrix has 2n rows and m = Y^a=i I^M columns, so = (0, . . . , 0) and 1 = (1, . . . , 1) 
denote row vectors of appropriate lengths. We introduce n new variables yi, ■ ■ ■ ,y n m order to 
encode the system /i = ■ • • = f n = in one polynomial with support in A, via the Cayley trick: 
4>{x,y) = yxhix) + ■■■ + y n fn(x). Note that %{A U . ..,A n ) = [Z 2n ,ZA\. 

In what follows when we refer to resultants or discriminants we will refer to the equations 
of the corresponding cycles, but we will omit the tildes in our notation. More explicitly, we will 
follow the convention in the article [5] by D 'Andrea and Sombra. In general, both definitions 
coincide, but this convention allows us to present cleaner formulas. For instance, when the family 
Ai, . . . ,A n+ i is essential, our notion of resultant equals the resultant in [T71 [25] raised to the 
index i(Ai, . . . ,A n+ i). In most examples these two lattices coincide, and so our resultant cycle 
equals the resultant variety and the associated resultant polynomial is irreducible. 

Remark 1. Assume A\ consists of a single point a and that {1} is the only essential subfamily 
of a given family A±, . . . , A n+ \. Let f\{x) = cx a . Then, for any choice of Laurent polynomials 
/*2 , . . . , fn+i with supports A2, ■ ■ ■ , A n+ i, it holds that (cf. J3, Proposition 2.2]) 



ResAi, 



,(/!,...,/„) = c MV(A 2 ,..,A n+1 )_ 



(3) 



With this convention, the following multiplicativity formula holds: 



1. [22] Let A^A'I, Ai, . . . , A n+1 be finite subsets of . 

"n+l 



Theorem 

/1, . . . , fn+i be polynomials with supports contained in A\ 
where f[ has support A[ and f" has support A'{. Then 



A 



i n with Ax = A[ + A'{. Let 
and assume that f\ = /{/{' 



ReSAi,...,i4 n+ i(/i ) • • • , fn+i) = Res A / iv .. iAn+1 (/{, . . . , f n+ i) ■ Res^ ... )j4n+1 (/", . . . , f n+x ). 

Cattani, Cueto, Dickenstein, Di Rocco and Sturmfels in [3] proved that the degree of the 
mixed discriminant A is a piecewise linear function in the Pliicker coordinates of a mixed Grass- 
manian. An explicit degree formula for plane curves is also presented in [31 Corollary 3.15]. In 
case A\,A2 consist of all the lattice points in their convex hulls, they are two dimensional and 
with the same normal fan, then the bidegree of A AljAa satisfies the following: bidegree of A Al)J 4 2 
in the coefficients of fi equals: 



Vol(Qi + Q2) — area{Qi) — perimeter (Qj), 



where i € {1,2}, i / j. where Qi = conv(Ai), i = 1,2, and Q\ + Q2 is their Minkowski sum. 
The area is normalized, so that a primitive triangle has area 1 and the perimeter of Qi is the 
cardinality of dQi n 1? . We will recover the general formula for this degree and present it in 
Corollary [5j 

Buse and Jouanolou consider in [2] the following equivalent definition of the mixed discrim- 
inant, in case where /1, . . . , f n are dense homogeneous polynomials in (xo, . . . , x n ) of degrees 
d\, . . . ,d n respectively, that is, their respective supports Ai = d^ are all the lattice points in 
the di-th dilate of the unit simplex a in R n . It is the non-zero polynomial in the coefficients of 
/1 , . . . , f n which equals 

for all i € {l,...,n}, where Jj is the maximal minor of the Jacobian matrix associated to 
/1, . . . , f n obtained by deleting the i-ih. We give a more symmetric and general formula in 
Corollary H] below. 

The multiplicativity property of the discriminant in the case of dense homogeneous polyno- 
mials was already known to Sylvester |26J and generalized by Buse and Jouanolou in [2J, where 
they proved that when in particular A\ = d\o = (d[ + d'{)a and f\ is equal to the product of 
two polynomials /{ • /{' with respective degrees d[,d'{, the following factorization holds: 

<^dia,...,dnO-{jli ■ ■ ■ 1 fn) = ^d' 1 o-,...,d n oA/l) • • • > fn) ' ^■d"a,...,d n ff\fl i • • • i fn) 

' R es d' 1 (T,d'{(T,...,dn<T(fl' fl > ■ • • ifn) ■ 

It is straightforward to see in general from the definition, that the vanishing of any of the 
polynomials A^ _ >An (f[, • • • , f n ), &a'(,...,aM"' ■■■Jn),oi Res A > iA ^_ An (f[, f", ...,/„) implies 
that 

^A' 1 +A'-[,...,A n (flfl, h, ■ ■ ■ , fn) = 0. 

It follows from |14] that when each support configuration Ai is full, the Newton polytope 
of the discriminant A-a' +A" ,A 2 ,—,An (/i/i'> hi ■ ■ ■ > fn) equals the Minkowski sum of the Newton 
polytopes of the discriminants A A / )A2) ...,A n (/i, h, ■ ■ ■ , fn) and &>A'{,Az,...,A n {fiJ2,---,fn) plus 
two times the Newton polytope of the resultant Res A / a" A 2 ... A„(f{i f'xi fii ■ ■ ■ ■> fn)- So, a first 
guess would be that the factorization into the three factors in ([5]) above holds for general supports. 
We will see in Corollary [6] that indeed other factors may occur, which we describe explicitly. 

This behaviour already occurs in the univariate case: 

Example 1. Let A' x = {0,ii, . . . ,i m ,di},A" = {0, j\, . . . , j/,^} be the support sets of f[ = 
o-o + a^x 11 + • • • + a,i m x lm + a^x^ , /{' = 60 + b^x^ 1 + • • • + bj t x^ 1 + bd 2 x d2 respectively. Then 

A(f[X) = A(f[).A(tt)-R(f[,n 2 .E, 

where E = a l Q~ m ° ft^ 1- " 10 a d 1-4m ~ m b£~ Jl ~ mi , with mo := mm{ii,ji} and mi := min{di — 
im,d2 — ji}- On the other hand, in the full case i\ = j% = 1, i m = d\ — l,ji = di — 1, thus E = 1 
because its exponents are equal to zero. 



3 A general formula 



The aim of this section is to present a formula which relates the mixed discriminant with the 
resultant of the given polynomials and their toric Jacobian, whose definition we recall. 

Definition 1. Let fi(xi, . . . , x n ), . . . , f n (xi, . . . , x n ) be n Laurent polynomials in n variables. 
The associated toric Jacobian jj equals x\ ■ ■ ■ x n times the determinant of the Jacobian matrix 
of f, or equivalently, the determinant of the matrix: 



Xi 



dx\ 



dfn 

. dxi 



dfi 

r -"n n 
OX„ 



dfn 



dx r 



Note that the Newton polytope of jj is contained in the sum of the Newton polytopes of 

Jit • • • t in- 
ks, we remarked before, we will mainly deal in this chapter with the case n = 2. Also, to 
avoid excessive notations and make the main results cleaner, we assume below that A\,Ai are 
two finite lattice configurations whose convex hulls satisfy 

dim(Qi) = dim(Q 2 ) = 2. 

Let /i,/2 be polynomials with respective supports A\,Ai: 



fi(x) 



a£Ai 



1,2, 



where x = (x\, x-z). We denote by £ the set of primitive inner normals n G (Z, 2 )* of the edges of 
A\ + A%. We call A^ the face of Ai where the inner product with n is minimized. We call this 
minumum value v T - '. We also denote by fj 1 the subsum of terms in /j with exponents in this face 



/*M 



E 



Ci.n>% 



1,2, 



which is rj- homogeneous of degree v%. Up to multiplying fi by a monomial (that is, after 
translation of Ai) we can assume without loss of generality that is? ^ 0. Now, A^ is either a 
vertex of Ai (but not of both A±,A2 since two vertices do not give a Minkowski sum edge), or its 
convex hull is an edge of Ai (with inner normal 77), which we denote by e^. Note that if the face 
of A\ + Ai associated to n is a vertex, both polynomials f^ are monomials and their resultant 
locus has codimension two. 

We denote by Hi(rj) (i = 1, 2) the integer defined by the following difference: 



Hi{rj) = min{(?7, m),m G Ai - A]} - z/' 



(6) 



and by 

fi(n) = min{/ii(77),/x 2 (r?)}, (7) 

the minimum of these two integers. Note that by our assumption that dim(Qj) = 2, we have 
that fj,(rj) > 1. 

Without loss of generality, we can translate the support sets A\,A\ to the origin and consider 
the line L v containing them. The residue (cycle) Res^ A^ifi-ifo.) ls considered as before, with 
respect to the lattice LTlZ 2 . 

Remark 2. As in Remark {J\ if f% is a monomial, the resultant equals the coefficient of f% 
raised to the normalized length t{e^) of the edge e\ of A2 (that is, the number of integer points 
in the edge, minus 1). If n is an inner normal of edges A\ and A\, then the resultant equals the 
irreducible resultant raised to the index of 7LA\ + 7LA\ in L'' fl Z 2 . In particular, the exponent 
/j,(rj) = 1 if at least one of the configurations is full. 

The following is our main result. 

Theorem 2. Let /i,/ 2 be generic Laurent polynomials with respective supports A\,A2- Then, 

R.esA 1 ,A 2 ,Ai+A 2 (/i,/2,'// ) = ^Ai,An(fl, f2) ■ E, 
where the factor E equals the finite product: 

E = X\Re SA n Al (flf 2 )^). 

r?e£ 

Proof. Let X be the projective toric variety associated to A1+A2. This compact variety consists 
of an open dense set Tx isomorphic to the torus (C*) 2 plus one toric divisor D v for each rj € E. 
The Laurent polynomials /1, / 2 , Jj define sections Li, L 2 , Lj of globally generated line bundles 
on X. The resultant ResA 1 ,A 2 ,Ai+A2(fii H-, Jj) vanishes if and only if L\,L2,Lj have a common 
zero on X, which could be at Tx or at any of the D^. 

There is an intersection point at Tx if and only if there is a common zero of /i,/2 and 
jj in the torus (C*) 2 . In this case, the discriminant A^4 1:J 4 2 (/i, / 2 ) would vanish. It follows 
that A J 4 1)J 4 2 (/i) f2) divides KesA 1 ,A 2 ,A 1 +A 2 {fi, Hi Jj)- (the indices [I? : 7LA\ + 2L4 2 ] and [Z 2 : 
ZAi + ZA 2 + Z(Ai + A 2 )} are equal). 

If instead there is a common zero at some D v , this translates into the fact that fi,f$ an d 
(Jj) 71 = Jj v (with obvious definition) have a common solution. But as f? are 77-homogeneous, 
they satisfy the weighted Euler equalities: 

mxi 7 ^ + mx2 7 ^ = v?fi, 2 = 1,2, (8) 

ax\ 0x2 

from which we deduce that Jj v lies in the ideal ^(/i,/ 2 ) and so, the three polynomials will 
vanish exactly when there is a nontrivial common zero of /-j* and f% ■ This implies that all facet 
resultants Res A v A v(f^, ff) divide Res A 1 ,A 2 ,A 1 +A 2 (h, f 2, Jj)- 



Now, we wish to see that the resultant Res^ A v Ui->f2) ra i se d to the power fj,(jj) occurs as 
a factor. The following argument would be better written in terms of the multihomogeneous 
polynomials in the Cox coordinates of X which represent L\,L 2 ,Lj [4J. Fix a primitive inner 
normal direction n € S of A\ + A 2 , let t be a new variable and define the following polynomials 

F t (t,x) = Y, c i>a t^-^x a , i = 1,2, (9) 



so that 

and we can write 



Fi(l,x) = fi(x), F i (0,x) = fV(x), z = 1,2, 



f?(x) = FV(t,x)-t^G i (t,x), » = 1,2, (10) 

where the polynomials G{ are defined by these equalities. The polynomials F\,F2, Jp define the 
sections L\,L 2 ,Lj. For each t, we deduce from the bilinearity of the determinant, that there 
exists a polynomial H(t, x) such that the toric Jacobian can be written as Jp = jT n +t ,J '^H(t, x). 
But, as we remarked, jj v lies in the ideal /(/j 7 ,/^), and using the equalities (fTOl) . we can write 
jj = Hi(t,x) + t^H 2 (t,x), with Hi € L(F 1 ,F 2 ). Note that if for instance r/i / 0, then the 
power of x\ in each monomial of Fi can be obtained from the power of t and the power of x 2 , 
that is, we could use t and x 2 as "variables" instead. We will denote by Res the resultant 
defined over X [3]. Therefore, 

ReSA^A^+A.ihJ^JJ) = Re4 uA2Al+ A 2 (Fi,F2,t^H 2 ). 

Now, it follows from Theorem Q] that 

is a factor of ResA 1 ,A 2 ,A 1 +A 2 (fi, f2, Jj)- Indeed, no positive power of t divides H 2 for generic 
coefficients. Considering all possible r] £ £ we get the desired factorization. □ 

Theorem [2] and the proof will be extended to the general n-variate setting in a forthcoming 
paper [?]. We only state here the following general version without proof. Recall that a lattice 
polytope P of dimension n in W 1 is said to be smooth if at each every vertex there are n 
concurrent facets and their primitive inner normal directions form a basis of Z n . In particular, 
integer dilates of the unit simplex or the unit (hyper)cube are smooth. 

Theorem 3. Let P C W 1 be a smooth lattice polytope of dimension n. Let Aj = (diP) n Z n , 
i = l,...,n, di,...,d n G Z>o> a- n d /l>---)/n polynomials with these supports, respectively. 
Then, we have the following factorization 

ReSA 1 ,..., J 4 n ,Ai+-+A n (/l,- ■■,fn,Jf) = &Ai,...,A n (fl, ■ • • , In) ' E, 
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where the factor E equals the finite product: 

E=l[Res A v !A v(flfZ). 

Note that all the exponents in E equal 1 and all the lattice indices equal 1. 

When the given lattice configurations Ai are the lattice points did of the dj-th dilate of the 
standard simplex a in W 1 , (that is, in the homogeneous case studied in [2j), formula (j4]) gives 
for any n in our notation: 

^■d 1 a,...,d n a(fl, ■ ■ ■ ,fn) ■ ^^(di(T) e i,...,(dnO-) e i (/l% ■ ■ ■ > fn)l 

where eo,... ,e n are the canonical basis vectors (or eo = — e\ — ■ ■ ■ — e n , if we consider the 
corresponding dehomogenized polynomials, by setting xq = 1). Note that Theorem [3] gives the 
following more symmetric formula: 

Corollary 4. With the previous notation, it holds: 

R es d 1 a,...,d n cr,(d 1 -\ h<2„)<r(/l> • • • > fn, Jf ) = 

n 

&diff,...,dnff{h> • ■ ■ s/n) ■ ||ReS( dl(T )< !ii ... i ( rfnCr )e i (/f ! ,. . . ,/^*). 

i=0 

It is straightforward to deduce from this expression the degree of the homogeneous mixed 
discriminant, obtained independently in [U|21|2T]. Similar formulas can be obtained, for instance, 
in the multihomogeneous case. 

We recall the following definition from [3] . If v is a vertex of Ai , we define its mixed multiplicity 

as 

mm AlA2 (v):=MV(Q 1 ,Q 2 )-MV(C i ,Q 3 ), {i,j} = {1,2}, (11) 

where Cj = conv(Ai — {v}). 

Let £' C E be the set of inner normals of A\ + A2 that cut out, or define, edges e^ in both 
Qi,Q2- The factorization formula in Theorem [2] can be written as follows, and allows us to 
recover the bidegree formulas for planar mixed discriminants in [3j. 

Corollary 5. Let A\,A<i be two lattice configurations of dimension 2 in the plane, and let 
/ij/2 be polynomials with these respective supports. Then, the resultant of/1,/2 and their toric 
Jacobian, namely ResA 1 , J 4 2 ,Ai+A 2 (/i) fii Jj), factors as follows: 

Aw/1,/2) • n c r Ai ' Aa(v) • n toB^u?,® 1 **- (12) 

v vertex of Ai or A 2 tjgE' 
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The bidegree (<5i,<52) of the mixed discriminant A^ 1) ^ 2 (/x, f 2 ) in the coefficients of fi and f 2 , 
respectively, is then given by the following: 

Vol(Q 1 ) + 2-MV(Q 1 ,Q 2 )-Y, i ( e ])-^)- Y, rnrn AuA2 (v), (13) 

r/£E' v vertex of (Ai) 

where i G {1,2}, i ^ j. 

Proof. To prove equality (|12p . we need to show by Theorem [2] that the factor 

r;e£ 

equals the product 

v vertex of A\ or A2 r)£.Y,' 

When r\ € X', i.e. r/ is a common inner normal to edges of both Qi, we get the same factor on 
both terms, since that our quantity /i(r]) equals the index min{it(ei(r/), A±), u(e2(r]), A 2 )}, in the 
notation of [3]. 

Assume then that r/ is only an inner normal to Q 2 . So, A\ is a vertex v, f% = ex" is a 
monomial (with coefficient c) and f% is a polynomial whose support equals the edge e\ of A 2 
orthogonal to r/. In this case, Res^ A v (fii /2O = ^ n ' by Remark [TJ 

For such a vertex u, denote by £(v) the set of those 77' ^ £' for which v + e\ is an edge of 
Q\ + Q2- Note that it follows from the proof of [31 Prop. 3. 13] (cf in particular Figure 1 there), 
that there exist non negative integers //(?/) such that 

mm 0) = Y^ t(e\ ) ■ fJ>'(v')- 
ij'e £(v) 

Indeed, /j(t/) = //(?/). 

To compute the bidegree, we use the multilinearity of the mixed volume with respect to 
Minkowski sum. Observe that the toric Jacobian has bidegree (1, 1) in the coefficients of fi, f 2 , 
from which we get that the bidegree of the resultant ResA 1 ,A 2 ,A 1 +A 2 {fi> /2, Jj) is equal to 

{2MV{A 1 ,A 2 ) + Yo\{Q 2 ), 2MV(Ai,A 2 )+Vol(Qi)). (14) 

Substracting the degree of the other factors and taking into account that the bidegree of the 
resultant Res^ A v {fiif 2 ) equals (^(e^),^^)), we deduce the formula (fT3l) . as desired. □ 
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4 The multiplicativity of the mixed discriminant 

This section studies the factorization of the discriminant when one of the polynomials factors. We 
make the hypothesis that /{, /{', / 2 have fixed support sets, and A^A", A 2 C Z 2 . So fi = f[- /{' 
has support in the Minkowski sum A\ := A\ + A"; in fact, its support is generically equal to A\. 
We will denote by /(/(r/) (resp. n"(rj)) the integer defined in ©, with A\ replaced by A\ (resp. 

A'{). 

Corollary 6. Assume A'^A'l and A2 are full planar configurations of dimension 2. Let /{, /{', $2 
be generic polynomials with these supports and let fi = f[- /{'• Then, 

A Al Mh »/a) = A^ifih) ■ A A ,, A2 (f",f 2 ) ■ Res A>lA2 (fifl f 2 f ■ E, 
where E equals the following product: 

II Res WiA3 ((/{)V 2 y' (?7) - M?7) ■Res K , M ((/nV 2 V' W - M?7) . (15) 

r?gS 

Proof. By Theorem [21 we get that 

a ft t \ R Ai,A2,Ai+M(fl,h,JJ) 

Al ' Aa(/l ' /2) " WRA^iflf,)^ ' (16) 

r,£S 

and similarly for A^/ ^ (/{ , / 2 ) and A^//^ (/{', / 2 ). Let us write the numerator of (fTBT) as follows: 

- R A' 1 +A' 1 ',A 2 ,A , 1 +A'/+A 2 (/l/l > /2i J f{f[',f 2 )' 
where JL„ , = fiJT„ t + fi'JT, t ■ Let us apply Theorem [Tito re- write it as follows: 

^A' 1 ,A 2 ,A' 1 +A' 1 '+A 2 (/l> /2> J f[f{',h) R A'{,A2,Ai+A'{+A 2 (fli h, J f[f['J 2 ) = 

= ^A'^Aa.Ai+A'Z+Aal/l, /2, f"Jf[,f 2 ) ^A' 1 ',A 2 ,A' 1 +A' 1 '+A 2 (./l / , /2, /l^/{',/ 2 )) 

because the resultant of {/ii, /i 2 + g/ii, . . . } equals the resultant of {hi, /i 2 , . . . }, for any choice 
of polynomials h±,h 2 ,g (with suitable supports). We employ again Theorem [1] to finalize the 
numerator as follows: 

-Ra' 1 ,A 2 ,A' 1 +A 2 (/i ; hi J f[,f 2 ) ' R A';,A 2 ,A"+A 2 (fl, /2, J f['J 2 ) ' ^ ,A" ,A 2 (^1 ' -^ ' ^) ■ 
For the denominator of (|16p . we use again Theorem [TJ to write: 

n ^-^(/r^r'^ • n R Ar,Ai(fPJ%y" {v) = 
U R A's + ArMf ,r inf v 2r {v) -E, 

r?gS 
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because the products 

II R A? ,AiU[\ttY' {ri) = II RA> r ,Ai(fr,f v 2r" M = h 

since fi^f^ (resp. fi V ,f%) are both monomials. To conclude the proof, simply assemble the 
above equations. □ 

As a consequence, we have deg A A2 A(/i, $2) = 

= deg^ A(/{, f 2 ) + deg^ A2 A(/f, / 2 ) + 2 • deg A , )A ,, jA2 i?(/{, f», f 2 ) - deg(E). 

When all the configurations are full and with the same normal fan, all the exponents h(v) = 
//(??) = /i"(f?) = 1. Therefore, E = 1 and no extra factor occurs. 

We define h'i (77) , h'{ (77) as in ^. Indeed, we now fix 77 and will simply write h'i-, ///, fx\, Hi- It 
happens that only one of the factors associated to 77 can occur in E with non zero coefficient. 
More explicitly, we have the following corollary, whose proof is straightforward. 

Corollary 7. With the notations of Corollary® for any 77 € £ it holds that: 

• If Hi = fl", then n' = jjl" = ji and there is no factor in E "coming from 77". 

• If Hi 7^ h'{> assume wlog that Hi = Mi < Mi- There are three subcases: 

~ /jfM2 < Mi j again there is no factor in E "coming from 77". 

— If Hi = Hi < M2 < Mi? ^ en ^ e resultant ReS( A /^ A '7((/{) 7? , /^ ) does no£ occur, but 
Res( A //\ t7jA '7((/{ / ) r ', /g ) /ias nonzero exponent (this resultant could just be the coefficient 
of a vertex raised to the mixed multiplicity). 

— If Hi = Hi < Mi < M2> ^e situation is just the opposite than in the previous case. 

5 Conclusion and future work 

The intent of this book chapter was to present our main results relating the mixed discriminant 
with the sparse resultant of two bivariate Laurent polynomials with fixed support and their toric 
Jacobian. On our way, we deduced a general multiplicativity formula for the mixed discriminant 
when one polynomial factors as / = /' • /". This formula occurred as a consequence of our main 
result, Theorem [2j and generalized known formulas in the homogeneous case to the sparse setting. 
Furthermore, we obtained a new proof of the bidegree formula for planar mixed discriminants, 
which appeared in [3]. 

The generalization of our formulas to any number of variables will allow us to extend our 
applications and to develop effective computational techniques for sparse discriminants based 
on well tuned software for the computation of resultants. 
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